Towards the identification of a quantum critical line 
in the (p, B) phase diagram of CeCoIn 5 
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The low-temperature thermal expansion of CeCoIns single crystals measured parallel and perpen- 
dicular to magnetic fields B oriented along the c axis yields the volume thermal-expansion coefficient 
/3. Considerable deviations of j3(T) from Fermi- liquid behavior occur already within the supercon- 
ducting region of the (B,T) phase diagram and become maximal at the upper critical field B" 2 - 
However, /3(T) and the Griineisen parameter F are incompatible with a quantum critical point 
(QCP) at £?J? 2 , but allow for a QCP shielded by superconductivity and extending to negative pres- 
sures for B < Bc2- Together with literature data we construct a tentative (p, B, T) phase diagram 
of CeCoIns suggesting a quantum critical line in the (p, B) plane. 

PACS numbers: 74.70.Tx, 73.43.Nq, 65.40.De 
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Many Ce-bascd hcavy-fcrmion compounds exhibit un- 
conventional superconductivity on the verge of long- 
range magnetic order, i.e. in the vicinity of magnetic 
quantum critical points (QCPs) [l|, ■ The presence of 
these magnetic instabilities is manifested by profound de- 
viations from Fermi-liquid (FL) behavior [3(. CeCoIns 
is a typical representative of these systems with pro- 
nounced non-Fermi-liquid (NFL) behavior at tempera- 
tures above its superconducting transition at T c = 2.3 K 
0] • With increasing magnetic field B the deviations from 
FL behavior become maximal at the upper critical field 
B c2 = B c2 (T -> 0) « 5T for B \\ c suggesting that NFL 
behavior might originate from a nearby field-induced an- 
tifcrromagnctic (AF) QCP at B® 2 - Indeed, the low- 
temperature dependence of the specific-heat coefficient 
C/T at B° 2 follows the Hertz-Millis-Moriya (HMM) the- 
ory 0-0] of a two-dimensional spin-density-wave insta- 
bility, as it might be expected for the layered, tetragonal 
crystal structure of CeCoIns 0, Q . The linear thermal- 
expansion coefficient along the c axis a c [jj], the resis- 
tivity, and the thermal conductivity fioj ]. however, re- 
veal at T < 0.3 K deviations from HMM behavior that 
were attributed to a dimensional crossover for T — > Q . 
Hall-effect measurements indicate a QCP at the second- 
order phase boundary of the recently discovered Q phase 
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12j , possibly indicating that critical fluctuations of 
the Q phase produce NFL behavior. In contrast, resis- 
tivity measureme nts p oint to a QCP at negative hydro- 
static pressures p 1 1 31 ] or emerging from a pressure- and 
field-dependent critical line that is masked by supercon- 
ductivity (l4j |. 

To search for the origin of NFL behavior of CeCoIns in 
the (p, B) plane, we have exploited a characteristic QCP 
feature: the accumulation of entropy at finite tempera- 
tures due to the instability of the ground state [l5| . It can 
be accessed from the derivatives dS/dp, or dS/dB, the 
volume thermal-expansion coefficient /3 = —V~ x ■ dS/dp 



(V = molar volume) and the magnetization dM/dT = 
dS/dB, respectively. Measurements of M point to quan- 
tum criticality directly at or below B® 2 (lij . In the vortex 
phase, however, M cannot be used to determine dS/dB 
due to irreversible flux motion. The thermal expansion, 
on the other hand, is less sensitive to pinning effects be- 
cause of the small nearly pressure-independent entropy 
of the vortex lattice. 

To study dS/dp, we performed thermal-expansion 
measurements in magnetic fields of up to 5 = 14 T 
with B || c between 40 mK and 4K. As an important 
extension to previous measurements [^, [T3| our capaci- 
tive dilatometer allows for thermal expansion and mag- 
netostriction measurements both parallel and perpendic- 
ular to the applied magnetic field with a resolution of 
Al = 10 -3 A. CeCoIns single crystals were grown in In 
flux. The plate-like crystals have typical dimensions of 
l a ~ 3 mm and l c w 0.5 mm thus enabling a relative res- 
olution of Al/l w 10~ 9 . To obtain the Griineisen param- 
eter T at high fields, we measured the specific heat C{T) 
between 350 mK and 5K at B = 10 T and 14 T with a 
Physical Properties Measurement System from Quantum 
Design. 

In a Fermi liquid the thermal-expansion coefficient ap- 
proaches a linear temperature dependence for T — > 0, 
proportional to the specific heat. The proportionality 
constant is essentially the Griineisen parameter T — 
V/k s ■ P/C (k s — adiabatic compressibility) and re- 
flects the volume dependence of the Fermi energy. To 
expose deviations from FL behavior in CeCoIns, we 
plot the volume thermal-expansion coefficient, calculated 
from the measured linear thermal-expansion coefficients 
/? = 2a a + a c , as /3/T in Fig.QJa) and (b) for fields below 
and above B® 2 , respectively. In zero magnetic field, /3/T 
drops asymptotically towards zero for decreasing T, just 
as C/T, due to the presence of the superconducting gap 
in the energy spectrum of the quasiparticles. At higher 
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fields, still in the field range of the second-order super- 
conducting phase transition at B C 2, /3/T grows and starts 
to diverge for T — > 0. To our knowledge, this is the first 
time that such a behavior is observed within the vortex 
phase of a superconductor, where, usually, the remanent 
Fermi liquid in the vortex cores leads to a finite, constant 
/3/T. This NFL behavior corresponds to that observed 
in specific-heat measurements where the entropy conser- 
vation requires a continuously increasing C/T ex — log? 1 
of the normal conducting background for T — > 0, Hf . 

(3/T exhibits the strongest low-temperatures upturn 
close to B® 2 [Fig. 0Ja)]. Further enhanced fields weaken 
the divergence and lead, in a growing temperature range 
at T < Tfl, to the restoration of the Fermi liquid, i.e. 
/3/T = const [Fig. 1(b)]. This clear feature does not 
support an interpretation Q in terms of a dimensional 
crossover. We emphasize that at low temperatures /3/T 
remains positive in the entire measured field range, ex- 
cept for the peaked anomaly of the first-order transition 
[Fig. QTc)]. This is an important result, because the ab- 
sence of a sign change demonstrates that S has no max- 
imum as a function of p at B® 2 but grows further under 
negative, hydrostatic pressure or, equivalently, with in- 
creasing unit-cell volume. These results speak against a 
QCP in the vicinity of B° 2 (p = 0). 

A direct proof for the existence of a pressure-tuned 
QCP is the Griineisen parameter T. If a system is domi- 
nated by a single characteristic energy scale E*(V), Y re- 
flects its volume dependence: T = — d\nE*/dln V. Since 
at a QCP E* vanishes, Y is expected to diverge when the 
QCP is approached [18J. We estimate Y by assuming k s 
to be roughly equal to the isothermal compressibility of 
CeCoIn 5 k t = 1.31 • 10" 2 GPa" 1 and with specific- 
heat data from [1,0, 21 1 and our measurements. The re- 
sulting T(T) curves, shown in Fig.[5Ja), exhibit high r(T) 
values up to 80, typical for heavy-fermion systems [22} . 
that illustrate their high sensitivity to volume changes. 
The T(T) curve at B = 5.1 T w B c ° 2 is the highest of all 
fields. Here, Y initially shows a rise with decreasing T, 
but levels off below w 260 mK and drops again. Although 
the strongly enhanced Y points to a pressure-tuned QCP, 
the lack of a low-temperature divergence provides clear 
evidence against a QCP at ambient pressure in the entire 
field range up to 14 T. 

To clarify a possible link between superconductivity 
and quantum criticality in CcCoIns , we compare the vol- 
ume dependence of T c with the measured Griineisen pa- 
rameter. If superconductivity provides the dominant en- 
ergy scale E* oc T c , Y is determined by — din T c /d In V — 
(T c k s ) _1 • dT c /dp at T « T c . The hydrostatic pres- 
sure dependence dT c /dp = 2dT c /da a + 2dT c /da c is cal- 
culated from the uniaxial pressure(tTi) dependences using 
the Ehrenfest relation dT c /doi = T C V • Ac^/AC (i = a, 
c). Here, Actj and AC are the discontinuities of the linear 
thermal-expansion coefficients and the specific heat at T c , 
respectively The resulting dT c / dp data are displayed as 
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FIG. 1. (Color online) The volume thermal-expansion coeffi- 
cient of CeCoIns divided by temperature /3/T in magnetic 
fields parallel to the c axis below (a) and above B® 2 (b). 
Note that the "overshooting" at T c is due to the fact that 
different samples with slightly different B C 2(T) curves were 
used to measure a a and ot c . The data for B = 4.85 T are a 
combination of measurements at B = 4.8 T and 4.9 T. The 
dashed lines represent the non-Fermi-liquid behavior for the 
HMM model, (c) /3/T as a function of B at T = 0.4 K ob- 
tained from magnetostriction measurements (continuous line) 
and from temperature-dependent measurements at constant 
B (circles). Dashed line is a guide to the eye. 



a function of B in Fig.[5Jb). They are in good agreement 
with experiments under hydrostatic pressure Q . dT c / dp 
shows a sign change at B + rj 4.1 ±0.2 T. At higher fields, 
when the transition becomes first-order, dT c /dp remains 
negative up to as demonstrated by the negative vol- 
ume discontinuities at T c [Fig.QJa)], which is expected to 
produce a likewise negative low-temperature Griineisen 
parameter in the field range B + < B < B® 2 ■ Experimen- 
tally, however Y(T) exhibits large positive values in the 
entire field range for T — > 0, apart from the small field 
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FIG. 2. (Color online) (a) The Griineisen parameter V versus 
T at different constant fields parallel to the c axis, (b) The 
hydrostatic pressure dependences dT c /dp as a function of B. 
Lines are guides to the eye. 



internal of the first-order transition near B C 2 [Fig. QJc)]. 
Therefore critical fluctuations of superconductivity or the 
first-order transition of the Q phase cannot be responsi- 
ble for the NFL behavior. 

The field dependence of the FL recovery [Fig. [TJb)] 
can be used to estimate a possible field-induced origin 
of the NFL behavior in CeCoIns at p = 0. To look for 
such a behavior with a decrease of Tfl toward a QCP, 
we construct from our data the (B, T) phase diagram 
in Fig- EJa). Within the experimental error, Tfl agrees 
with the onset of the T 2 dependence of the resistivity 
P EH HI or of constant C/T values 0). Other char- 
acteristic features that are related to the FL recovery 
such as the minimum of the differential Hall coefficient 
25| , coincide likewise with Tp l ■ A striking feature of the 
phase diagram is that Tfl does not vanish at B^ 2 - For 
B < 10 T, Tfl extrapolates linearly to a T = critical 
field of B c = 4.1 ± 0.1 T, which is equal to B + where, 
according to dT c /dp = 0, T c exhibits a maximum as a 
function of p. This suggests that the (hidden) AF QCP 
at p = is expected to occur at the field where T c (p) of 
the superconducting dome is maximal. Furthermore, the 
evolution of the coherence maximum of the magnetore- 
sistance p(B) with temperature, T(p max (B)), [Fig. [3ja)] 
[HI extrapolates to roughly the same B c for T —> 0. 
Therefore a QCP at B c has been postulated |25| . In the 
T = plane with two control parameters p and B, a line 
of QCPs is highly plausible. This is supported by the 
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FIG. 3. (Color online) (a) The (B, T) phase diagram of 
CeCoIns at ambient pressure constructed from our measure- 
ments of the superconducting transition (•), the onset of FL 
behavior at Tfl (a), and the change in the critical behav- 
ior at T cr (▼). In addition, literature data are shown for the 
maximum of the magnetoresistance T(p m ax(B)) (o [24[) and 
Tfl determined by Hall effect (□ [25|) and resistivity mea- 
surements (o [2(1). (b) Tentative (p, B, T) phase diagram 
of CeCoIns. The data from (a) have been extended by pub- 
lished measurements under hydrostatic 0, [27|- [29| and under 
negative, chemical pressure on CeCoIns - -r Cd^ |30l - [33 ] . To 
account for the different sample qualities and varying T c def- 
initions the measurements have been scaled to match in the 
regions of overlap. B+ (•), and the onset of the NFL behavior 
(■ [3) are projected onto the T = plane. 



thermal-expansion measurements above T c [Fig. QJa)] or 
above Tfl [Fig.QJb)] where all j3/T curves merge into the 
same critical T _1 dependence. Therefore, the NFL be- 
havior at different fields is likely to have a common origin. 
The NFL behavior might be explained by a single QCP 
at negative critical pressure p c and B B+ or, alterna- 
tively by an extended quantum critical phase boundary. 
Indeed, high-pressure experiments indicate that the NFL 
behavior emanates from ap- and B-dependent line, which 
is shielded by superconductivity [l4j . 
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To check whether a quantum critical line is compati- 
ble with the p dependence of T c and the NFL behavior 
from our measurements, we extend the phase diagram to 
finite pressures with literature data measured under hy- 
drostatic pressure or on Cd-doped CeCoIn 5 [Fig. EJJb)]. 
Although Cd doping does not generate a noteworthy 
enlargement of the unit cell, pressure experiments on 
CeCoIn5_. J ;Cd 2: show that a nominal Cd concentration 
of 5% correponds to a negative chemical pressure of 
— 0.7GPa |33|. Doping with Cd suppresses superconduc- 
tivity and leads to long-range AF order. At B — 0, a 
linear extrapolation of the Neel temperature Tn to pos- 
itive pressures reveals that, without intervening super- 
conductivity, AF order would disappear at a pressure 
where T c becomes maximal as indeed shown above for 
B = B c k, B + . This is a common feature of many un- 
conventional superconductors close to AF order and is 
generally taken as evidence for superconductivity medi- 
ated by magnetic fluctuations 0, HI- As shown in the 
T = plane of Fig. Et^b), with increasing B, the T c max- 
imum follows approximately the origin of NFL behavior 
to lower pressures until it reaches at p = a field that is 
close to B + (or B c ). At higher fields, especially at B® 2 , 
the phase line continues to negative pressures, again in 
accordance to our findings. In contrast to the p and B 
dependence of the NFL behavior, the onset of the discon- 
tinuous superconducting transition and Q phase B n (p) 
follows the upper critical field B® 2 (p) and is still present 
when, beyond p > p c , NFL behavior can no longer be 
observed. Although, at ambient pressure, Bo and B + 
are close to each other, their qualitatively different pres- 
sure dependence gives clear evidence that the Q phase 
can be ruled out as source for the NFL behavior. The 
close relationship between the T c maximum, the origin 
of NFL behavior, and the endpoint of AF order, on the 
other hand, suggest that the quantum critical line cor- 
responds to the AF phase boundary. In this case, the 
available data can be described by the tentative phase 
diagram proposed in Fig. [3jb) . 

In summary, we performed thermal-expansion mea- 
surements at different magnetic fields to study the pres- 
sure and field dependence of the non- Fermi-liquid behav- 
ior in CcCoIns which extends at ambient pressure over a 
wide field range. Although the deviations from Fermi- 
liquid behavior reach a maximum at the onset of su- 
perconductivity at B® 2 , the Griineisen parameter clearly 
demonstrates that at p = no QCP exists up to a field 
of 14 T. Due to the finite distance of the QCP from B® 2 , 
the critical behavior changes at low temperatures toward 
Fermi-liquid recovery. A combination of our results with 
literature data allows to construct a tentative (p,B,T) 



phase diagram, in which the NFL behavior originates 
from a quantum critical line of the onset of antiferro- 
magnetic order. This line, emanating from a point of 
positive pressure on the B = axis, is hidden by the 
superconducting dome. It may extend to negative pres- 
sures at magnetic fields larger than B® 2 as suggested by 
recent experiments on Cd-doped CeCoIns. 
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